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Abstract 

We outline a generalization of the Grobner fan of a homogeneous ideal with maximal cells 
parametrizing truncated Grobner bases. This "truncated" Grobner fan is usually much smaller 
than the full Grobner fan and offers the natural framework for conversion between truncated 
Grobner bases. The generic Grobner walk generalizes naturally to this setting by using the 
Buchberger algorithm with truncation on facets. 

We specialize to the setting of lattice ideals. Here facets along the generic walk are given by 
unique (facet) binomials. This along with the representation of binomials as integer vectors give 
an especially simple version of the generic Grobner walk. 

Computational experience with the special Aardal-Lenstra integer programming knapsack 
problems is reported. The algorithms of this paper are implemented in the software package 
GLATWALK, which is available for download at http : //home . imf . au . dk/ni el s /GLATWALK 

1 Introduction 

The generic Grobner walk Q is a version of the classical Grobner walk algorithm for Grobner basis 
conversion in the Grobner fan of an ideal in a polynomial ring. In the generic walk explicit rational 
vectors in the Grobner fan are replaced by computations with infinitesimal numbers which can be 
handled formally. This leads to an algorithm where input consists only of a source Grobner basis, a 
source term order and a target term order. 

Truncation of homogeneous ideals have proved very valuable in algebraic computations related 
to integer programming (see for example |7 1). In general there are a lot fewer truncated initial ideals 
than initial ideals. Similary to initial ideals, truncated initial ideals may be parametrized by the 
maximal cells in a complete polyhedral fan. We introduce this fan, which is easily constructed from 
the usual Grobner fan by inserting the truncation operator at the appropriate places. The truncated 
Grobner fan is in general much smaller than the full Grobner fan and forms the polyhedral setting 
for a truncated version of the generic Grobner walk. We prove that the truncated Grobner fan is 
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regular along the lines of |6|. This leads to a "truncated" state polytope with vertices enumerating 
the different reduced truncated Grobner bases. 

In the setting of lattice ideals we give a rather detailed version of the generic Grobner walk. 
Algebraic computations with lattice ideals can be greatly simplified representing (saturated) bino- 
mials as integer vectors. This along with the fact that the generic walk only traverses facets lead to 
several simplifications. We report on computational experience in computing saturations of lattice 
ideals and truncated test sets related to the integer programming problems posed in Q. Our exper- 
iments show that the the generic walk in the truncated Grobner fan consists of significantly fewer 
steps, whereas the walk in the full Grobner fan does not compare well computing directly with the 
Buchberger algorithm. 

I am grateful to B. Sturmfels for inspiring conversation and for greatly simplifying my original 
approach to truncated Grobner fans. Thanks are also due to K. Fukuda, A. N. Jensen and R. Thomas 
for slowly making me grasp the joys of algebra and polyhedral geometry. 

2 Preliminaries 

We let R = k[x\, . . . ,x n ] denote the ring of polynomials over a field k. We will view R as the semi- 
group ring fc[N n ]. 

2.1 Grading on R 

Given n elements a\, . . . ,a n of an abelian group (A, +) we let Sa denote the semigroup N«i H h 

Na n C A. For v = (vi, . . . , v n ) G N" we put deg(x v ) = v ■ a = V\a\ H h v n a n . These data give a 

natural A -grading on R by defining 

R s = span k {x v | deg(x v ) = s} 

for s G Sa- A non-zero element / G R is called homogeneous of degree deg(/) = s if / G R s . Given 
an ideal J in R, we let J s = J n R s - Recall that an ideal / is homogeneous ideal if J = Q)J S and that 
this is equivalent to J being generated by homogeneous elements. We call R positively graded if 
Ro = k. This is equivalent to dim k R s < °° for every s G Sa- 

2.2 Truncating subsets 

A subset Q C 5a is called truncating if s,t G Q. whenever s + 1 G CL for s,t G Sa- Our standard 
example of a truncating subset is 

Q.b = {x G Sa I b — x G Sa} 
for b G Sa- To a truncating subset Q. C Sa we associate the monomial ideal (cf. Remark 12.8 in [6|) 

Ma = (x v | deg(v) Q) CR. 
Given a homogeneous ideal J C R we let 

se£i 
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Lemma 2.1 Let I and J be homogeneous ideals in R. Then 



la = Jq. if and only if I+M& = J + M a . 



Proof. If / is a homogeneous ideal, then I + M& is a homogeneous ideal and 




This proves the lemma. 



□ 



2.3 Initial ideals 

Let -< denote a total multiplicative ordering on monomials in R (we do not require that the monomial 
1 is minimal). If ft) G W 1 , we let ~< m denote the multiplicative ordering defined by x 11 ~< m x v if 
co- u < co vorco u = co v and x u -< x v . For f £ R, we let supp(/) denote the set of monomials 
occuring with non-zero coefficient in /. We let in^(/) denote the maximal (initial) term in supp(/) 
with respect to -<. Similarly we let in ffl (/) denote the sum of terms a v x v in supp(/) with co ■ v 
maximal. For a subset G C R we let in ffl (G) and in^(G) denote the ideals (m a (f) \ f G G\ {0}) and 
(m~>(f) | / G G) respectively. These ideals are homogeneous if / is homogeneous. A Grobner basis 
for / over -< is a finite set G := {/i, . . . ,/ r } C / such that 



The Grobner basis G is called minimal if none of f\, . . . ,f r can be left out and reduced if in^(/,) 
does not divide any of the terms in /,■ for i ^ j and i,j= 1 , . . . , r. The reduced Grobner basis of an 
ideal is unique and consists of homogeneous elements if the ideal is homogeneous. A homogeneous 
ideal J in R always has a reduced Grobner basis over -< if R is positively graded, since dim^7 s < °° 
for s G Sa- We record the following simple but crucial result (|6|, Proposition 1.8) with a complete 
proof. 

Proposition 2.2 Let I C R be any ideal and CO G E n . Then 

in < (m a {I))=m <m {I). 

Proof. Clearly in^ ffl (7) C in^(in ffl (/)). The ideal in m {I) is homogeneous in the grading given by 

CO. So we may decompose an element / G in^/) as / = H \- fx,> where fx. G ina,^) is 

homogeneous of ft)-weight Xj for j = 1, . . .,t. Now in^(/) = in^(/^.) for some j and we may write 



Notice that the multiplicativity of -< is not used in the proof of Proposition l2.2l The lifting from 
in ffl (7) to / in the proof is a key element in the Grobner walk algorithm. 



in^(/) =in^(G) = (mMi),---^Mr))- 



f x . = a x in ra (/i ) + ••• + a r va. a (f r ) 
for suitable f\ , . . . , f r G /, where a\ , . . . , a r are homogeneous elements. Therefore 

hj =in(B(«l/l H VClrfr) 

and in^(/) = va. <m {a\f\ H Ya r f r ). This shows that in Ko) (7) D m^(m a (I)). 



□ 
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2.4 Truncated Grobner bases 

Let 7 be a homogeneous ideal in R and Q. a truncating subset. A finite subset G C 7 is called an 
Q-Grobner basis for 7 over -< if 

in^(7) fl = in^(G) Q . 

If the coefficients of the initial terms in^ (g) are 1 for g G G and in^ (g) does not divide any of the 
terms in g' for g ^ g' G G, then G is called a reduced Q-Grobner basis. Reduced Q-Grobner bases 
are unique. 

Proposition 2.3 Let G be the reduced Grobner basis for J over -<. Then G consists of homogeneous 
elements and Gq = {g G G | deg(g) G £1} is the reduced Q-Grobner basis for J over -< 

Proof. If g is an element of the reduced Grobner basis of 7 and g = g\-\ h g r is written as 

a sum of homogeneous elements, then in^(g) = 'm^(gj) for some j = l,...,r. Therefore g has 
to be homogeneous. The monomial ideal in^(7) is spanned as a vector space by {in^(/) | / G 
7 S , for some s G S}. This shows that in^(7)n is the fc-span of in^(/) for / G 7^- Suppose that 
/ G Jq. Since G is a Grobner basis for 7 we may find g G G, such that in^(g) divides in^(/). This 
shows that deg(g) G £2, since Q. is a truncating subset. Therefore in^(/) G (in^(g) | g G Gq)q. □ 

Corollary 2.4 7/" G^ is reduced Q.-Grdbner basis for I over -< a then {m m {g) \ g G Gq} is the 
reduced Q.-Grobner basis for in^/) over -< 

Proof. Let G be the reduced Grobner basis for I over Then we know from Proposition l2.2l that 
{m a (g) | g G G} is the reduced Grobner basis of in ft) (7). Now Proposition l2.3l gives that 

{^M I g e G,<teg(infl,(g)) G a} = {in ffl (g) | g G Gn} 

is the reduced ^-Grobner basis of /. □ 

Proposition 12.31 reveals that the ^-truncated Grobner basis can be obtained from the reduced 
Grobner basis by picking out the elements with degree in £1. Truncated Grobner bases can be com- 
puted from a homogeneous generating set using Buchbergers algorithm discarding 5-polynomials 
with degree outside Q.. This follows from the fact that the division algorithm preserves the degree of 
a homogeneous polynomial. Very often only Grobner bases up to a certain degree are needed. 

3 The truncated Grobner fan 

In this section we assume that R is positively graded. Analogously to ([6 1, Proposition 2.3) we define 

Ca[fi>] = {v G R n | in v (/) Q = m a (r) a } 
for a homogeneous ideal I C R. We call the closure of Cn[a>] in R w a truncated Grobner cone. 
Theorem 3.1 The collection 

&£l(I) = {Qi[w] I CO G M"} 
of truncated Grobner cones form a complete fan in M. n . 



5 



Proof. A monomial ideal M satisfies m u (J)+M = in M (7 + M) for every ideal J C R and u G R". Now 
Lemma l2~Tl shows that in v (/)a = in a (I)ci if and only if m v (I+Mo) = in ffl (/+M^). This proves that 
#b(7) is the usual Grobner fan of the homogeneous ideal I+M&. Now the conclusion follows from 
(0, Proposition 2.4). □ 

The truncated Grobner fan is available from the usual Grobner fan by eliminating inequalities 
given by polynomials of degree outside £1. We give an example illustrating this. 

Example 3.2 Consider the (toric) ideal 

Ia = (a 2 c — b 2 e, a 2 d — be 2 , ce — bd) C k[a, b : c, J, e\ . 

This ideal is homogeneous in the grading given by the columns of 

/I 1 1 1 1\ 
A= 1 2 1 . 
\0 1 2 1/ 

For example, the degree of the variable c is (1,2, 1). The Grobner fan &{Ia) ofI& is the normal fan 
of an octagon in R 5 (cf. Example 1.1 in It is pictured in R 2 below with reduced Grobner bases 
labeling the maximal cells. 
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Putting ii = {vG Sji | v ■ ( 1 , 1 , 1 ) < 6} we get the truncated fan ^q(/a) with the reduced Q-Grobner 
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bases labeling the maximal cells. 




4 Truncated state polytopes 

The truncated Grobner fan of / is the usual Grobner fan of I + M&. From this it follows that the 
truncated Grobner fan of a homogeneous ideal / is the normal fan of a natural Minkowski summand 
in a state polytope for /. Emphasizing the simple Lemma lOl below. we briefly sketch a proof of this 
along the lines of |6|. Let / be a homogeneous ideal in R and a G A. 

Lemma 4.1 Let -< and -<' be two total multiplicative orderings on monomials in R. Let x Vl , . . .,x Vs 
denote the monomials in m^(I) a and x 111 , . . . ,x Us the monomials in in x > (I) a . Then after permuting 
u\ , . . . , u s and v\,...,v s we may assume that x 1 ' 1 -< x Vl -< ■ ■ ■ -< x Vs and 

x 11 ' -< x Vl 
x" 2 -< x V2 

x lls ^x Vs . 

Proof. We may find a vector space basis f\ , . . . , f s of I a , such that in^ (f\) = x Vl , . . . , in^ (f s ) = x Vs 
and x Vl -< x V2 -<■■■-< x Vs . Now put f[ = f\. Move on to f 2 . If in^(f 2 ) = Ain^/(/{) for A G k, 
put f 2 = f 2 - Xf[. Then clearly inu</(y£) -< in^(/ 2 )- In general if in^(/J) -< in^(//) for j < m and 

in^/(/ m ) 6W = span A .{in^/(/ 1 / ), . . tnen /« = fm - Mf[ hn-\f n -\ satisfies 

m ^'(f!n) & W for suitable X\,. . . , A m _i G fc. Furthermore in^/(//J -< in^ (/„,). In this way we get the 
monomials x l<1 = in^/ (/{),..., x" s = in^ of in^/ (7) a written up in the desired way. □ 

Let 

Ein_<(/) a = £ v 

x v em^(I) a 

where a G A and -< is a multiplicative total ordering. 
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Corollary 4.2 Let -<\ and -<2 be two total multiplicative orderings. IfL in^, (I) a = E in^ 2 (I) a , then 
in Kl (/) a = in^ 2 (/) a . 

Proof. This is an easy consequence of Lemma |4~T1 □ 

Definition 4.3 A state poly tope in degree a for I is defined as 

State a (7) = conv{E in^(/) a |-< total multiplicative ordering} . 
For a finite subset S C Sa we let 

State 5 (/) = £State fl (/). 

Corollary 4.4 For CO eR n we have 

face ffl (State fl (/)) = State^in^/)). 

Proof. First assume that in ffl (7) is a monomial ideal and face ffl (State a (/)) is a vertex {E in^(/) a } for 
some multiplicative monomial ordering -< by picking a generic ft). Apply Lemma l4~Tl to in^ ffi (I) a and 

in^(/) a . In this setting we then have u\ =v\,...,u s = v s since ft) - (vi H hv s ) > ft) - {u\ H 1-"*). 

Therefore {E in^(/) a } = {E in£o(/) a } when ft) is generic. General ft) are reduced to generic ft) as in 
the last part of the proof of Lemma 2.6 in 1 6 1 . □ 

Now let G be a universal Grobner basis for / consisting of homogeneous elements. Put S = 
{deg(g) \g EG}. Given a truncating subset Q C 5a we let 

State n (7) = StatennsW 

denote a truncated state polytope. Notice that State^/) is a Minkowski summand in a state polytope 
States (7) of I. Using Corollaries 14.21 and 14.41 the same arguments as in the last part of the proof of 
Theorem 2.5 in show that the normal fan of States (I) is J^h(/). It follows that J^hCO is the 
normal fan of a Minkowksi summand of a state polytope and that J^h(7) is a coarsening of the usual 
Grobner fan. 



5 Walking in the truncated Grobner fan 

The Grobner walk can be carried out in the truncated Grobner fan converting one truncated Grobner 
basis to another. We sketch the appropriate generalization of Proposition 3.2 in [3|. The term orders 
of (Q, Proposition 3.2) are represented by weight vectors below. If a weight vector ft) represents 
the term order -<, then we refer to ^ as ft) modified by 77 . 

Proposition 5.1 Let I be a homogeneous ideal in R and Q. a truncating subset. Let Qj[ft>i] and 
Cfl[o>2] be maximal cells in the truncated Grobner fan ^ofj) of 1. Suppose that G is the reduced 
Q.-Grdbner basis for I over (0\. If CO G Co[ff>i] HQifft^], then 

(i) The reduced Q.-Gr6bner basis for in a (I) over C0\ is G m = {in co (^) | g e G}. 
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(ii) IfH is the reduced Q.-Grobner basis for in a (I) over C0i, then 

{f-f G \feH} 

is a minimal Q.-Grdbner basis for I over CO2 modified by CO. 

( Hi) The reduced Q.-Grdbner basis for I over Oh modified by CO coincides with the reduced Q.- 
Grdbner basis for I over C0i- 

Proof. The items (i) and (iii) follow as in Proposition 3.2 of Q taking Proposition ^ . 31 and Corollary 
12 .41 into account. For the proof of (iii) observe that if H' = {f\, . . . ,f s } is the reduced Grobner basis 
of in,y(7) over CO2, then G" := {f\ — /p , . . . ,f s — ff } is a minimal Grobner basis for / over C0i 
modified by co by Proposition 3.2(ii) in Q. Here G' is the reduced Grobner basis for / over C0\. 
Notice that G" consists of homogeneous elements and that 

G'l 1 = {geG"\dcg( g )ea} 

is a minimal Q-Grobner basis for / over CO2 modified by CO. If deg(/ ; — fp') G £1, then deg(/ ; ) G £1. 
In this case ff' = ff. Using Proposition 12 . 31 this finishes the proof of (ii). □ 

Now the generic Grobner walk (Q, §4) carries over verbatim to the truncated setting using the 
Buchberger algorithm with truncation in step (iv). 

6 Lattice ideals 

In the rest of this paper we will remain exclusively in the setting of lattice ideals. Recall the decom- 
position of an integral vector v G Z" into v = v + — v~, where v + , v~ G N" are vectors with disjoint 
support. For u, v G N n we let u < v denote the partial order given by v — u G N n . For a subset SS C Z" 
we associate the ideal 

I@=(x v+ -x v ~ \ ve-^)cR. 

In the case where SS = Jzf is a lattice we call the lattice ideal associated to Jzf . Recall that lattice 
ideals are saturated i.e. if / G I^f is divisible by a variable xu then f/xi G I^f. This means that we 
apply the homogeneous Buchberger algorithm with sat-reduction as explained in 0. 
Define 

bin(w) =x w+ -x w ~ 

for w G Z n . The (saturated) 5-polynomial of bin(i<) and bin(v) is then given by bin(« — v). Similarly 
if v + < u + we may reduce bin(w) by bin(v) giving bin(w — v). We have silently assumed that the 
initial term of bin(w) is x w for the term ordering in question. 

Usually a generating set for Jzf as an abelian group is given. Computing the lattice ideal 
Il£^>Im can be done using that 

I^ = lg§- (Xl---Xn)°°, 
where / : f°° denotes the ideal given by 

{reR\ rf m el, form>0} 

for an ideal I C R and an element / G R ([6|, Lemma 12.2). If SB contains a positive vector, then 
I 38 = 1% (0' Lemma 12.4). If Jzf fl N" = {0}, 1% may be computed from lag using Grobner basis 
computations for different reverse lexicographic term orderings ([6|, Lemma 12.1). 
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6.1 The generic Grobner walk for lattice ideals 

We now specialize the generic Grobner walk to the setting of lattice ideals representing binomials 
by integer vectors as bin(v). In our implementation of the generic walk we walk from -< CJ to -< C2 , 
where -< is the reverse lexicographic order given by x\ -< • • • -< x n and c\ ,c 2 are integer vectors. 

In the algorithm outlined below we walk between two arbitrary multiplicative orderings -<i, -<2 
and move inside Grobner cones given by minimal Grobner bases (cf. Q, Proposition 2.3). Autore- 
duction is replaced by a simplified lifting step. 

In the notation of (Q, §2.3) we have 

<5^(bin(w)) = w 

assuming that x w -< x w+ . The facet preorder -< is now given on binomials bin(a) and bin(v) as in 
(□,§4, (3)) by 

bin(w) -<bin(v) •<=>- TuV -<i Tvu' , 

where T is a matrix defining the target term order -<2- We get as in ([3|, §4) that bin(w) -< bin(v) 
and bin(v) -< bin(w) imply that u is a multiple of v. If G := {bin(vi), . . . ,bin(v r )} is a Grobner basis 
then vi, . . . , v r lie in a common half space and the facet preorder induces a total ordering on G. In 
particular one gets that the facet in the generic Grobner walk is given by a unique facet binomial. 
This means that Grobner basis computations on facets proceed as in ([4J, Algorithm 3.1). To give 
some more details we introduce the notation 

mon(w) = x w+ 

for w £ Z". Suppose that G above is a minimal Grobner basis and that bin(vi ) is minimal in the facet 
preorder. Then step (c) of ( iPfl. Algorithm 3.1) is to compute a Grobner basis of the ideal 

(bin(-vi),mon(v 2 ),. . .,mon(v r )}. (1) 

Working with minimal Grobner bases it may happen that mon(— vi) is divisible by a monomial 
mon(vj) , j = 2, . . . , r. In this case the reduced Grobner basis of © is 

{mon(vi ),..., mon(v r )} 

and mon(vi) lifts to bin(w), where bin(— w) is the reduction of bin(— vi) modulo 

{bin(v 2 ),...,bin(v r )}. 

In this way the usual autoreduction of the Grobner walk is built into the lifting. Notice that mon(v 7 ) 
lifts to bin(vy) for 7 = 2,..., r. These observations account for step (ii) in facetJbuchberger below. 
On the other hand, if mon(— vi) is not divisible by any of the monomials, then bin(vi) is a (real) 
facet binomial of a facet in the Grobner cone corresponding to the reduced Grobner basis. In this 
case we end up with a minimal Grobner basis 

{bin(— vi),mon(vt>i), . . .,mon(w s )} 

of the ideal in CO). This lifts to the minimal Grobner basis 

{bin(-vi),bin(wi), . . . ,bin(wj)}. 
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The details of the algorithm are given below. The variable facet Jist contains a list of binomials 
ordered in ascending order according to the facet preorder (these are potential facet binomials). The 
variable G contains the current minimal Grobner basis. The procedure initialize_facet_list initializes 
G and facet Jist given B. The procedure insert inserts a given binomial into G and updates facet Jist. 

Notice that we do not really compute the S-polynomials in (iv.b) below. We ODtimize the algo- 
rithm by replacing the 5-polynomial 5(bin(w) , mon(v) ) with the initial term x^ v ^ of the saturated 
5-polynomial sat (S (bin (w) , bin(v) ) . 

Algorithm 6.1 (Generic Grobner walk for lattice ideals) 

INPUT: Integer vectors c\,C2- Integer vectors B = {vi , . . . , v r } such that {bin(vi ) , . . . , bin(v r ) } is a 
minimal Grobner basis for I ^ over -<\. 

OUTPUT: Integer vectors G = {w\ , . . . , w s } such that {bin(wi ) , . . . , bin(w ? ) } is a minimal Grobner 
basis over -<2- 

( i) initialize_facet_list, 

(ii) while (facet Jist 7^ 0) do 

(a) facet-bin :— first element in facet Jist 

(b) facetJbuchberger, 

facetJbuchberger: 

(i) Delete f acet Join from G and put bin := — facet Join; 

(ii) ifw + < bin + for some w G G 

reduce bin by G; 
insertf -bin); 
return; 

(Hi) Spairs := 0; 

(iv) for v in G do 

(a) if(bin + Av+ = 0) 

continue; 

(b) Spairs := Spairs U {v — bin} 

(v) Delete v G G ifbin + < v + ; 

(vi) while (Spairs j^&) do 

(a) Select s in Spairs and put Spairs := Spairs \ {s}. 
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(b) Reduce s by bin; 

(c) if (v + < s + ) for some v G G 

continue; 

(d) Delete vEG ifs + < v+. 

(e) inserts 

(f) S pairs : = S pairs U {s — bin} 
(vii) insert(bin); 

Truncation blends in easily with Algorithm 16. 11 Suppose that Q. denotes the truncating subset. 
First binomials in B with degrees outside Q. are discarded. With every addition of an 5-binomial in 
step iv(b) of facetJbuchberger, a test for degree membership of the truncating subset Q. is done. If 
the test fails for the 5-binomial it is not added to Spairs. 

Example 6.2 We give a very simple example illustrating Algorithm \6. 1\ Consider the ideal 

1= (x-t 2 ,y-t 3 ) C k[t,x,y]. 

Clearly G := {x — t 2 , y — t 3 } is a Grobner basis for I over the weight vector o = ( — 1,0,0), where t 
corresponds to (1,0,0) etc. We wish to walk to the weight vector % = (1,0,0) breaking ties with the 
reverse lexicographic order given by t < x < y. Let -< denote the corresponding facet preorder. Then 
x~t 2 <y — t 3 and we begin by "computing" a Grobner basis for (t 2 —x,y) giving G= {t 2 — x, y — t 3 } 
after lifting. In the following step the facet binomial is y — t 3 , which gets replaced by the reduction 
y — tx in step (ii) o/ facetJbuchberger. This accounts for the next facet binomial. We then compute 
a Grobner basis of (tx — y,t 2 ) giving {tx — v,f 2 ,?y, y 2 }. This lifts to {tx — y,f 2 — x, ty — x 2 ,y 2 — x 3 }, 
which is the reduced Grobner basis for I over (1,0, 0), since ty — x 2 and y 2 — x 3 are not candidates 
for facet binomials as the vectors (1,-2,1) and (0,-3,2) both are outside C< CT ,< T (cf. §4 of 13]). 

7 Computational experience 

In Q a collection of integer knapsacks are constructed related to the classical Frobenius problem of 
finding the largest number, which is not a sum of given relatively prime natural numbers. Feasibility 
for these knapsacks turn out to be very hard for traditional branch and bound software like CPLEX, 
but easy for lattice reduction methods as shown in Q]. 

In these knapsacks are equipped with a feasible right hand side and a specific cost vector 
c. These examples form the point of departure in this section, where we specifically document 
performance for computing (truncated) test sets using the package GLATWALK 1 . It turns out that test 
sets in the feasibility case is by far the hardest computations. Test sets with respect to the cost vector 
in [2J finish in negligible timings (< 0.05 seconds) using both the generic walk and the Buchberger 
algorithm with truncation. 

Each of the examples are of the form: maximize cx, where 



Ax = b : 

'home . imf . au . dk/niels/ GLATWALK 



(*) 
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iGN" and A is a 1 x n-matrix (a\ ■ ■ -a n ). The cost vector c and the matrices A may be found in ||2l . 
The first step is finding a feasible solution to (*). As in Q this results in the knapsack: minimize 
t subject to Ax + 1 = b, where t E N and x E N'\ This leads to the problem of finding a (truncated) 
Grobner basis of 

( Xl -t a \...,x n -t a ») (**) 

with respect to the vector x = (1,0, ... ,0), where t is the "first" variable. We may compute this 
Grobner basis directly using the Buchberger algorithm or walk from the vector o = ( — 1,0, ... ,0). 
The performance of the functions walk and gbasis of GLATWALK for computing a full Grobner basis 
of (**) over x is reported in Q. The second step is the computation of the toric ideal Ia (associated 
with the integer matrix A) and its Grobner basis over the vector — c. The function saturate of 
GLATWALK performs the saturation necessary in computing lattice ideals. Below 2 we have computed 
the ideals Ia using saturate after LLL-reducing Ker(A) with the function LLL. In many of the 
examples, LLL-reduction offers great savings in the computation of the saturation. The third column 
shows the timing of gbasis in computing a full Grobner basis over — c for Ia- The fourth column is 
the timing of walk in walking from —e\ to — c. The fifth and sixth columns show sizes of the full 
and truncated reduced Grobner bases of Ia over —c. 



EXAMPLE 


saturate 


gbasis 


walk 


size 


tr size 


cuwwl 


0.1 


3.1 


1.8 


2618 


7 


cuww2 


0.0 


0.3 


1.3 


898 


16 


cuww3 


0.1 


0.4 


10.1 


963 


16 


cuww4 


0.0 


3.3 


59.6 


3143 


5 


cuww5 


0.0 


0.0 


102.8 


267 


32 


probl 


0.0 


0.0 


2.7 


180 


75 


prob2 


0.0 


0.0 


0.5 


280 


45 


prob3 


0.5 


0.0 


2.5 


163 


94 


prob4 


0.2 


0.1 


122.3 


475 


83 


prob5 


0.6 


0.0 


0.0 


68 


56 


prob6 


0.4 


9.2 


39.5 


4541 


94 


prob7 


0.4 


1.8 


72.8 


2036 


79 


prob8 


0.9 


0.0 


2.4 


227 


103 


prob9 


0.0 


0.0 


0.5 


108 


108 


problO 


1.4 


0.1 


517.5 


536 


119 



Both walk and gbasis finish in negligible timings (< 0.05 seconds) in computing the truncated 
Grobner bases. However in computing the full Grobner bases in the above table, walk does not 

2 All timings are in seconds. The computations were carried out on an ACER notebook 1.6 GHz Pentium mobile with 
1MB L2 cache. 
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compare well with gbasis. Typically to compute a target Grobner basis with less than 1,000 bino- 
mials, the generic walk traverses cones associated with reduced Grobner bases of more than 20,000 
binomials along a straight line intersecting many cones in the Grobner fan. 

Most of the examples above indicate that the truncated Grobner fan is much smaller than the full 
Grobner fan. The straight line path in the truncated Grobner fan traverses significantly fewer cones. 
It is open for further research exactly when the walk is a substantial improvement (as in many of the 
feasibility examples reported in Q). Perhaps a combination of direct Grobner basis computations 
for suitably chosen (easier) weight vectors tending to the target vector followed by a walk to the 
target order may lead to improvements. 
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